Interaction of a TeV Scale Black Hole with the Quark-Gluon Plasma at LHC 



Andrew Chamblin 1 , Fred Cooper 2 , and Gouranga C. Nayak 2 
1 Department of Physics, Queen Mary, University of London, Mile End Road, London El 4NS 
2 T-8, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA 

(February 1, 2008) 



m 

O 
O 

(N 
C 

a 
oo 



> 
On 

m 

o 
m 
O , 

^ : 

a: 
i 

Q . 
^ ■ 



X 



If the fundamental Planck scale is near a TeV, then parton collisions with high enough center-of-mass 
energy should produce black holes. The production rate for such black holes has been extensively 
studied for the case of a proton-proton collision at y/s = 14 TeV and for a lead-lead collision at sfs 
= 5.5 TeV at LHC. As the parton energy density is much higher at lead-lead collisions than in pp 
collisions at LHC, one natural question is whether the produced black holes will be able to absorb 
the partons formed in the lead-lead collisions and eventually 'eat' the quark-gluon plasma formed 
at LHC. In this paper, we make a quantitative analysis of this possibility and find that since the 
energy density of partons formed in lead-lead collisions at LHC is about 500 GeV/fm 3 , the rate of 
absorption for one of these black holes is much smaller than the rate of evaporation. Hence, we 
argue that black holes formed in such collisions will decay very quickly, and will not absorb very 
many nearby partons. More precisely, we show that for the black hole mass to increase via parton 
absorption at the LHC the typical energy density of quarks and gluons should be of the order of 
10 10 GeV/fm 3 . As LHC will not be able to produce such a high energy density partonic system, the 
black hole will not be able to absorb a sufficient number of nearby partons before it decays. The 
typical life time of the black hole formed at LHC is found to be a small fraction of a fm/c. 



I. INTRODUCTION 



There is now a huge literature devoted to the possi- 
bility that the scale of quantum gravity could be as low 
as a TeV [1]. If this is true, then quantum gravity and 
perhaps even string theory effects will become relevant 
in collider experiments. For example, there are many 
discussions on graviton and radion production and pro- 
cesses related to them at LHC [2,3]. One of the most 
exciting aspects of this TeV scale gravity will be the pro- 
duction of black holes in particle accelerators. These 
'brane-world' black holes will be our first window into 
the extra dimensions of space predicted by string theory, 
and required by the brane-world scenarios that provide 
for a low energy Planck scale [1]. While the exact met- 
rics describing black holes in brane-world scenarios are 
unknown, considerable work on this issue is underway 
[4]. Furthermore, even without the exact metrics it is 
possible to make estimates based on crude information. 
In particular, it is well-understood that when the mass of 
the black hole is greater than the Planck scale, the gravi- 
tational field of the brane can be neglected; furthermore, 
as long as the size of the black hole is small compared to 
the characteristic length scales, then a brane-world black 
hole may be regarded, to very good approximation, as 
simply a higher-dimensional black hole in flat space. Us- 
ing these approximations, in a number of recent papers 
people have studied the production of microscopic black 
holes in proton-proton (pp) and lead-lead (PbPb) colli- 
sions at LHC and cosmic ray events [5-16], [17], [18-20]. 

The evolution of a black hole after its formation is gov- 
erned by two competing effects: 1) absorption of nearby 
particles by the black hole and 2) evaporation of the black 



hole. Many papers have been written for the case of pri- 
mordial black holes within the braneworld scenarios [21]. 
In this paper we will study the case of the black holes 
formed at LHC. In addition to any black holes formed at 
LHC in pp or PbPb collisions, there will be many more 
standard model particles formed, in particular a huge 
number of quarks and gluons. The energy density of the 
quarks and gluons formed in a pp collision is not very 
large and hence the absorption of partons by the black 
hole is very small. However, the energy density of the 
partons formed in the PbPb collisions at LHC (at yfs — 
5.5 TeV per nucleon) is large. This is because a huge 
number of partons are formed in PbPb collisions (equiv- 
alent to 208-208 pp collisions) in a small volume. There- 
fore one needs to determine if a black hole can absorb a 
sufficient amount of partons after its formation in PbPb 
collisions to eventually 'eat' the quark-gluon plasma. In 
this paper we will analyze the evolution of the black hole 
at LHC by incorporating the absorption of partons by the 
black hole and the decay of the black hole simultaneously. 

The paper is organized as follows: In section II we 
present the calculation of the black hole mass evolution 
at LHC. In section III we present a perturbative QCD cal- 
culation to determine the energy density of the partons 
formed in PbPb collisions at LHC. In section IV and V 
we present our main results and conclusions respectively. 



II. EVOLUTION OF A TEV SCALE BLACK 
HOLE AT LHC 

The absorption rate of the black hole is proportional to 
the surface area of the black hole and to the energy den- 
sity of the mass less gluons produced at LHC [22]. The 
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effective black hole radius r e // for capturing particles is 
given by [23]: 
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where M P is the TeV scale Planck mass, M BH is the 
black hole mass and d is the number of extra dimensions. 
The effective radius in various numbers of extra dimen- 



sions is given as shown: 
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The cross section for the particle absorption by the 
black hole is given by: 



a = nr 



eff 
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Using this absorption cross section the accretion rate of 
the black hole mass becomes [22]: 
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accr = Firr% ff e(T), 
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where e(r) is the energy density of the mass less gluons 
produced at LHC. As gluons are the dominant part of 
the total parton production at LHC (more than 80 per- 
cent) we will consider the absorption of the gluons in this 
paper for simplicity. The energy density e(r) of the glu- 
ons formed in PbPb collisions at LHC will be calculated 
in the next section by using pcrturbative QCD. For an 
upper limit on our estimates on the accretion rate we will 
take F = 1 in our calculation. This value might be lower 
indicating less accretion [22] . As we are interested in the 
maximum accretion we will use F = 1 in this paper. 

Now we consider the evaporation rate of the black hole 
at LHC. ft is known that the evaporation is much faster 
for smaller mass black holes [24] and hence the rate of 
evaporation of the TeV scale black hole at LHC is much 
higher. We use the recent calculations of black hole radi- 
ation in the context of TeV scale physics [21,24,25] which 
yields the evaporation rate: 



dM 



BH i 



2 o ( 



120 

I (d±3\(d+3)/(d+l) 



(d+l) d + 3 



where 



7680 4^(^1)1 rl ff 



4nr 



BH 



(5) 



(6) 



with tbh being the Schwarzschild radius of the black 
hole: 
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and the area in the extra dimension is given by: 
(1) 
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As the change in Stefan-Boltzmann constant is very 
weakly dependent on the number of extra dimension [24] 
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we have used the four dimensional value equals ^ in 
the above. In the equation (5) g e ff corresponds to the 
effective degrees of freedom of the standard model par- 
ticles into which the black hole evaporates. The effec- 
tive number of degrees of freedom is given by: g e ff = 
<7bos + | flterm- In this paper we will assume that 
a black hole radiates mainly to the mass less particles 
(electrons and positrons, muons and anti-muons, pho- 
tons, neutrions, gluons and up, down, strange quarks) 
for which g e ff = 46 with <? e ^=18 for spin s = 1 and 
28 for s = |. T s is the dimensionless grey body factor 
which we take the same value as in [23] in the geometric 
optics approximations. There has been recent progress 
on grey body factors in extra dimensions [26]. ^2 in the 
above equation is the sum of fermions and bosons. For 
thermodynamics to be a reasonable approximation one 
requires Mbh >> Mp [23], which will be true for the 
ratio being > 3. 

The net change of the black hole mass is therefore: 



dM BH (T) dM BH ( T ). 
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We get / 2 =9.45, 



./4=30 and /7=210. By using the relation between r e ff 
and the black hole mass M BB we obtain, from the above 
equation: 



dM BH (T) 
dr 



= F d (r)M^(r) - 



B d 



(11) 



This equation has to be solved numerically once we know 
the gluon energy density e(r) at LHC. In the above equa- 
tion 
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We now just need to determine the gluon energy density 
at LHC by using perturbative QCD methods. 
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III. GLUON ENERGY DENSITY AT LHC 

Now let us determine the energy density of the mass 
less gluons formed in PbPb collisions at LHC. We use per- 
turbative QCD methods to determine the gluon energy 
density produced at LHC. In the lowest order of pQCD 
calculations the cross section for gluon-minijet produc- 
tion in pp collisions is given by [27-29]: 

Ojet = J dp T dy idy 2 —r^Yl x 1 fi/A{xi,P%) 



ijkl 



X2 fj/A(x2,Pr) <rij->ki(s,i, 



(13) 



The jet cross section can be related to the total number 
of minijets formed in PbPb collisions at LHC: 
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jet 



= T(0) J dp T dy 1 dy 2 —p^Yl Xl ^/a{xi,Pt) 

ijkl 



x 2 fj/A{x2,Pr) &ij^ki{s,i,u), 



(14) 



where T(0) = 9A 2 /8irR A is the nuclear geometrical fac- 
tor for head-on PbPb collisions (for a nucleus with a 
sharp surface). Here Ra = 1.2 A 1 / 3 fm is the nuclear ra- 
dius, (iij^ki denotes the elementary pQCD parton cross 
sections, which are given by: 
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As gluons are the dominant part of the parton produc- 
tion at LHC we have considered the above two partonic 
type collisions. The kinematic relations are given by: 



s = X1X2S = \-p\ cosh 2 



2/2/1-2/2 



(17) 



The momentum rapidities y\ , y 2 and the momentum frac- 
tions x\, x 2 are related by, 

xi = p T {e Vl + e V2 )/ Vi, x 2 = pt {e~ yi + e~ V2 ) /<J~s. 

(18) 

We use the saturation argument to fix the minimum mo- 
mentum scale as po above which the pQCD is applicable 
to be 2 GeV at LHC [30]. The partonic structure func- 
tions fij/A used in Eq. (13) are GRV98 [31] for the proton 
structure function and EKS98 [32] for nuclear shadowing 
effects. 

The total transverse energy < E^ 1 > produced in nu- 
clear collisions is given by: 

< E 3 T et >= T(0) f dp T p T d yi dy 2 ^^J2 x i fi/A{xi,p 2 T ) 



ijkl 



X2 fj/A(X2,Pr) &ij^kl(s,i,u). 



(19) 



The above pQCD formula has the information in the mo- 
mentum space via px and y (the momentum rapidity) 
distribution. In order to obtain energy density we need 
to supply additional information about coordinate dis- 
tribution. For an expanding system the minijet number 
distribution eq. (14) can be related to the phase-space 
distribution function via [33] : 



et 
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= 9c J f{x,p) p^da^ 
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where gc = 16 is the product of spin and color degrees 
of freedom, 



da^ = ir R A rdr] (cosh r], 0, 0, sinh??), 



(21) 



and 



p^ = (pt cosh y,PT cos <j>,pr sin (j>,pT sinhy), (22) 

for an 1+1 dimension expanding system. Using the above 
relations we get at the initial time tq (=l/po) : 
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where we assume a transverse isotropy at the early stage. 
We take a boost non-invariant distribution function in 
the following. We assume a Gaussian coordinate rapid- 
ity 7] (=[±ln^]) distribution of the form [29] 
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where f(pr, y) is obtained from by using pQCD 

from eq. (14). Using the above form in eq. (23) we get: 
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which gives: 
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From the above equation we get: 
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which gives a boost non-invariant initial phase-space 
gluon distribution function: 
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Using the above boost non-invariant gluon distribution 
function the initial minijet energy density is given by: 

£(to,v)=9c fdr(p' J 'u IJ ,) 2 f(pT,v,y,To) = gc j d 2 p T p 2 T 

j dy cosh 2 {q - y)f(pr,V,y, T o) = 2tt J dp T pr 
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Before we calculate the initial gluonic energy density 
at LHC we would like to find out the values of the un- 
known parameter a 2 appearing in the above equation. 
We will fix this unknown parameter by equating the first 
Pt moment of the distribution function with the pQCD 
predicted Et distribution as given by eq. (19). We have 

J dp T j d V jyJ pT = 2gcn 2 R 2 A T J dp T J dy p 3 T f{p T ,y) 



= / dp T / dy 



dN 3Ct 

^ * dydpT 



y/n<r(p T )e' 72 (PT)/ i ' 



(30) 



We determine the value of a 2 by using the pQCD values 

dE^ et 

of dy J pT from eq. (19) in the left hand side of the above 
equations. Assuming a 2 independent of pr we find from 
the above equation a 2 = 0.28 at LHC. Using this values 
of the a 2 we find from eq. (29): e =517 GeV/fm 3 at 
LHC in the central rapidity region. 



IV. RESULTS AND DISCUSSIONS 

In this section we will present the results for the evo- 
lution of the black hole mass at LHC by taking both 
accretion and radiation into account. The accretion rate 
depends on the energy density of the gluons (see eq. (11)) 
which is a function of time. The energy density of the 
gluons at LHC decreases as a function of time. This is be- 
cause of the fast expansion of the system, which acquires 
more volume. In our numerical calculation we will as- 
sume a maximum energy density at all time which is the 
energy density obtained at initial collision, as described 
in the last section. Thus our results will overestimate 
the effect of accretion and thus underestimate the crit- 
ical energy density. However the expansion time is at 
a fm/c scale which is slow compared to the evaporation 
time scale at LHC. We solve the eq. (11) numerically 
with the gluon energy density obtained in the last sec- 
tion. The time evolution of a black hole of mass 3 TeV is 
shown in Fig.l. The Planck mass is choosen to be 1 TeV. 
Different curves are for different numbers of extra dimen- 
sions. Solid, dotted and dashed lines correspond to d=2, 



4 and 7 respectively It can be seen that the decay of a 
3 TeV mass black hole is almost instantaneous at LHC 
even if one takes into account the absorption of partons 
with energy density of 517 GeV/fm 3 . The black hole of 
mass 3 TeV (with Planck mass equals to 1 TeV) decays 
within a fraction of a fm/c. The partons in the QGP 
phase hadronize in a hadronization time scale which is 
more than a fm/c. Hence a 3 TeV mass black hole com- 
pletely evaporates before any of the hadrons are formed 
at LHC. 

Now let us study the evolution of a heavier black hole, 
say 5 TeV, at LHC. As stated earlier a heavier mass black 
hole will evaporate slower. This is because the temper- 
ature of the black hole is inversely proportional to its 
mass. The evolution of a 5 TeV mass black hole at LHC 
is shown in Fig. 2. We have used a Planck mass equals 
to 1 TeV. In comparison to Fig. 1 the decay time of the 5 
TeV black hole mass is slightly enhanced when compared 
to the case when the mass of the black hole is 3 TeV, but 
it is still instantaneous in comparison to hadronization 
scale which is more than a fm/c. This suggests that a 
typical black hole formed at LHC evaporates instanta- 
neously and can not eat sufficient matter. 

In Fig. 3 we have shown the results of the time evo- 
lution of the 5 TeV mass black hole but with Planck 
mass equals to 2 TeV. As can be seen from the figure 
the evaporation is much faster than the earlier two cases, 
considered in Fig. 1 and 2. This is because as the Planck 
mass is increased the Schwarzschild radius is decreased. 
Due to the decrease in the Schwarzschild radius the effec- 
tive cross section for absorption of matter by the black 
hole is decreased. Hence the evolution of the black hole 
mass decreases much faster as we increase the Planck 
mass which can be seen from Fig. 3. 

In the following we will present estimates of the mini- 
mum energy density of the matter we would need in order 
for a TeV scale black hole to swallow quark-gluon matter 
at LHC. In order for a black hole to increase its mass, 
one must have: 



dM 



BH 



dr 



> 0. 



(31) 



Hence from eq. (10) we get the minimum required energy 
density (e min ) of the partonic matter to be: 



76807rr 4 c/ ' 



(32) 



From the above equation we get for M P = 1 TeV: 

1) Case with the number of extra dimensions d=2: 

e min = 4.2 x 10 10 GeV/fm 3 for M BH = 3 TeV 
e min = 2.1 x 10 10 GeV/fm 3 for M BH = 5 TeV, 

2) Case with the number of extra dimensions d=4: 
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e min = 3.2 x 10 10 GeV/fm 3 for Mbh — 3 TeV 
e min = 2.2 x 10 10 GeV/fm 3 for M BH = 5 TeV, 

3) Case with the number of extra dimensions d=7: 

e m in = 7.5 x 10 10 GeV/fm 3 for M BH = 3 TeV 
e mi n = 5.8 x 10 10 GeV/fm 3 for Mbh = 5 TeV. 

It can be seen that the minimum energy density required 
for different values of Mbh and for different values of 
number of extra dimension is much higher than the en- 
ergy density obtained at LHC. Note that these numbers 
are obtained for Mp= 1 TeV. For higher values of Mp 
these numbers will be further increased. Creation of such 
astronomically high values of energy density is beyond 
the reach of any accelerator experiment. This means that 
black holes created in high energy hadronic colliders at 
LHC will not swallow the matter nearby. 

V. CONCLUSIONS 

If the fundamental Planck scale is near a TeV, then 
parton collisions with high enough center-of-mass en- 
ergy should produce black holes. The production rate 
for such black holes has been extensively studied for the 
case of a proton-proton collision at y/s — 14 TeV and for 
a lead-lead collision at yfs = 5.5 TeV at LHC. As the 
parton energy density is much higher at lead-lead col- 
lisions than in pp collisions at LHC, it is necessary to 
check if the produced black hole will be able to absorb 
the partons formed in the lead-lead collisions and even- 
tually consume much of the quark-gluon plasma formed 
at LHC. In this paper, we have made a quantitative anal- 
ysis of this possibility and found that as the energy den- 
sity of partons formed in lead-lead collisions at LHC is 
about 500 GeV/fm 3 , the rate of absorption of these black 
holes is much smaller than the rate of evaporation. Hence 
black holes formed in such collisions evaporate much too 
quickly to absorb the partons nearby. For the black hole 
mass to increase via parton absorption at the LHC the 
typical energy density of quarks and gluons should be of 
the order of 10 10 GeV/fm 3 . As LHC will not be able to 
produce such a high energy density partonic system, the 
black hole will not be able to absorb sufficient nearby par- 
tons before it decays. The typical life time of the black 
hole formed at LHC is found to be a small fraction of a 
fm/c. 
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Black Hole Evolution at LHC 

M P =1 TeV, M EH =3 TeV 




Black Hole Evolution at LHC 

M P =1 TeV, M EH =5 TeV 



FIG. 1: Evolution of black hole of mass equals to 
3 TeV at LHC. The Planck mass equals to 1 TeV. 



FIG. 2: Evolution of black hole of mass equals to 
5 TeV at LHC. The Planck mass equals to 1 TeV. 



Black Hole Evolution at LHC 

M P =2 TeV, M EH =5 TeV 
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FIG. 3: Evolution of black hole of mass equals to 
5 TeV at LHC. The Planck mass equals to 2 TeV. 
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